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Abstract 

We present an autonomous phase-plane describing the evolution of Friedmann- 
Robertson- Walker models containing a perfect fluid (with barotropic index 7) in Brans- 
Dicke gravity (with Brans-Dicke parameter u>). We find self-similar fixed points corre- 
sponding to Nariai's power-law solutions for spatially flat models and curvature-scaling 
solutions for curved models. At infinite values of the phase-plane variables we recover 
O'Hanlon and Tupper's vacuum solutions for spatially flat models and the Milne uni- 
verse for negative spatial curvature. We find conditions for the existence and stability 
of these critical points and describe the qualitative evolution in all regions of the (to, 7) 
parameter space for < 7 < 2 and uj > —3/2. We show that the condition for infla- 
tion in Brans-Dicke gravity is always stronger than the general relativistic condition, 
7 < 2/3. 



1 Introduction 



One of the simplest extensions to the Einstein-Hilbert action of general relativity is the 
introduction of a scalar field non- minimally coupled to the metric tensor The simplest 
scalar-tensor theory is that proposed by Brans and Dicke @ where the gravity theory contains 
only one dimensionless parameter, uj, and the effective gravitational constant is inversely 
proportional to the scalar field, $. This theory yields the correct Newtonian weak- field 
limit, but solar system measurements of post-Newtonian corrections require u > 500 ||. In 
the limit u — > 00 the field $ becomes fixed and we recover Einstein gravity. This has led to 
more general scalar-tensor gravity theories 0, |5|, § being considered with a self-interaction 
potential V(Q) [0, §] or a variable u(&) || in order to fix $ by the present day. 

The fact that such scalar-tensor gravity theories are a generic prediction of low-energy 
effective supergravity theories from string theory [ 10 1 or other higher-dimensional gravity 
theories pLl has led to considerable interest in cosmological models of the very early universe 
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derived from Brans-Dicke gravity. The time-dependent gravitational constant introduces a 
new degree of freedom into cosmological models which has led different authors to propose 
extended inflation [12], pre big bang [Vjl and gravity- driven inflation models [14] of the early 
universe, as well as modifications to conventional inflation [I5j or nucleosynthesis models jl6j . 

In this paper we will investigate the qualitative evolution of Friedmann- Robertson- Walker 
(FRW) models containing a barotropic perfect fluid in Brans-Dicke gravity. In Section |3] we 
review the asymptotic behaviour for spatially flat FRW models and the general solution 
which can be given in parametric form for barotropic fluids [[H], Analytical solutions 
for spatially curved models are only possible in a limited number of special cases such as 
for radiation |Tj| [2(J pl| [14] or stiff fluid |20], [22| and other methods are required to study 



the general behaviour. In Section ^ the field equations describing the cosmological evolution 
of the Brans-Dicke field and scale factor are reduced to two coupled autonomous first-order 
equations. A phase-plane analysis [23] is then used to solve qualitatively for the cosmological 
evolution [^, |25|. We present for the first time a compactified phase-plane which represents 
both late and early time behaviour for flat (k = 0) and curved (k ^ 0) FRW models. In 
section [5] the cosmological evolution is described qualitatively for all regions of parameter 
space (7, u), where the Brans-Dicke parameter uj > —3/2 and the barotropic index of the fluid 
< 7 < 2. In Section 6 we consider whether it is possible to construct cosmological solutions 
with infinite proper lifetimes in some parameter regimes. Our results are summarised in 
Section 7. 



2 Equations of motion 

We will consider here the simplest scalar-tensor theory as originally proposed by Brans and 
Dicke, where the gravity theory is characterised by a single dimensionless parameter ou, 
representing the strength of coupling between the scalar field and the metric. The equations 
of motion are derived by extremizing the action 



S 



167T 



d 4 x sf—g 



(2.i; 



Homogeneous and isotropic cosmologies can be described by the Friedmann-Robertson- 
Walker (FRW) metric, 



-dt 2 + a{tf 



dr' 



1 — kr 2 



+ r 2 (d6 2 + sin 2 



(2.2) 



The evolution equations for the scale factor, a(t), and the Brans-Dicke field, $(£), can then 
be written as 
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where a dot denotes differentiation with respect to cosmic time, t. Matter with density p 
and pressure P obeys the continuity equation 



P = -3-(p + P). (2.5) 
a 

These evolution equations possess a first-integral which corresponds to the generalised Fried- 
mann constraint for the Hubble expansion, H = a/a, 

i^ = ^-4 + -(-V-4- (2-6) 
3$ $ 6\$j a 2 V ; 

Considerable insight can be gained into the dynamical evolution of these models by 
working with the conformally rescaled Einstein metric [^, |27j 

g„ v = (G N <S>)g^ , (2.7) 

where Gat is Newton's constant which remains fixed in the Einstein frame. The rescaled 
scale factor becomes a = (Gjy^Y^a and the cosmic time 

t = J (G N <5>) 1/2 dt . (2.8) 

The constraint equation (|2.6j ) then can be written as the familiar Friedmann constraint 

H 2 = —^(p + p)--, (2.9) 

o CL 

where H is the Hubble expansion in the Einstein frame, p = p/(Gn&) 2 , and 



is the effective energy density of the Brans-Dicke field in the Einstein frame. This has an 
effective pressure P = p, but is non-minimally coupled to the barotropic fluid. Note however 
that from the Friedmann constraint equation (|2.9|) we can immediately deduce that only 
closed models (k > 0) have a turning point H = in the Einstein frame for 2uj + 3 > 0. 
Theories with 2u + 3 < have a negative energy density in the Einstein frame, implying 
that the Minkowski vacuum spacetime is unstable, and so we shall not consider such models 
here. 

An equation of state for matter, P(p), is needed in order to solve these equations of 
motion. We restrict our analysis to barotropic fluids where the pressure is simply proportional 
to the density, P = (7 — l)p, which includes the case of pressureless dust (7 = 1), radiation 
(7 = 4/3) or a false vacuum energy density (7 = 0). The continuity equation ( ^5|) can then 
be integrated to give the matter density as a function of the cosmo logical scale factor a(t) 

P = ^z (2.11) 

a 61 
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and this allows us to eliminate both p and P from the remaining equations of motion. 

General analytic solutions for scalar-tensor cosmologies exist for only a limited number 
of special cases. For example we can solve analytically for a general FRW models in vacuum 
(p = 0) ]28|, ^l|] or when the matter is radiation (P = p/3) |20], [H], [14| or a stiff fluid 
(P = p) |20], 22] when the system reduces to perfect fluids in the Einstein frame. In Brans- 
Dicke gravity the general analytic solution for other cases can only be given in parametric 
form for spatially flat models (k = 0) and we shall review these results briefly before going 
on to study the evolution of the more general FRW models. 



3 Spatially- flat solutions 

Brans and Dicke's original motivation was to produce a model where the effective gravita- 
tional constant decreases ($ increases) as the density of the universe decreases. In their 
original paper [Q] , Brans and Dicke gave a power-law solution for the evolution of a spatially 
flat FRW model filled with pressureless dust 

a oc , $ oc It-t*! 9 , (3.1) 

where p = (2u;+2)/(3u;-|-4), and q = 2/(3u;-|-4). This was generalised to arbitrary barotropic 
equation of state by Nariai who found ||29| 



2u;(2 - 7) + 2 . , 

P = ^7 —7— — 3.2 

F 3u;7(2-7)+4 v ; 

= 2(4-37) 3) 

H 3cu7(2- 7 )+4 v ; 

This yields the general relativistic solution, a oc t 2//37 with $ constant, in the limit uj — > 00. 

However Nariai's power-law solution is in fact only a particular solution. The general 
solution can be given analytically in terms of a rescaled time coordinate 



Then the scale factor and Brans-Dicke field evolve as [1171 IIS 



where 



a = a (T -T^) n *(T -T + ) n± , (3.5) 
$ = $ (T-T_) m T(T-T+) m± , (3.6) 



n± = — (3.7) 

3(l + u;(2-7)=FVl + 2u;/3) 



1 ± y/1 + 2a,/3 

m± = (3.8) 

1 + oo(2 - 7) t V 1 + 2^/3 



The solution where T + = T_ reduces to Nariai's power-law solution given by Eqs. (|3.1|-|3.3|), 
and this is clearly the attractor solution when T — > 00. However, substituting a solution 
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of the form given in Eq. (|3.1|) into the definition of the time coordinate T in Eq. ( |3.4|) we 
can see that t —>■ oo at a finite value of T if 3(7 — l)p > 1. Thus Nariai's solutions are not 
necessarily the attractor solutions at late cosmic times. This case needs to be studied more 
carefully to determine the attractor behaviour which we shall do in the phase-plane analysis 
that follows. 

As T approaches T + or T_ the general evolution (for T_ 7^ T + ) approaches the vacuum 
solutions (p = 0) first presented by O'Hanlon and Tupper P5| , which are also power-law 
solutions of the form given in Eq. (|3.1|) with 



ti± oj + 1 ± y/l + 2u/3 

P± ~ 1 + 3(7- l)n± ~ 4 + 3c; ' 

m+ —2 

Q± = 7 x — = ; • (3-9) 

l + 3(7-l)n± l±3,/l + 2cj/3 



Nariai's solution is sometimes described as the matter-dominated (T — > 00) solution as 
opposed to the vacuum-dominated (T — > T±) solution. The vacuum-dominated evolution 
has two solutions which have been termed fast or slow |17| depending on the early time 
behaviour of the scalar field <£>. The slow solution corresponds to {p~,q-) defined above 
where at early times the scalar field $ is increasing (and G — > 0) as cosmic time increases 
for uj > —4/3. The fast solution is defined by (p+, q + ), with $ decreasing (and G — > 00) as 
t increases for uj > —4/3. 

The Brans-Dicke effective action Eq. (|2.1|) reduces to the truncated string effective action 
in the absence of matter (p = 0) when the coupling constant u = —1 [[TTJ. For comparison 
with the string literature we will define, 

= -ln$, a = lna. (3.10) 

It can be shown that the spatially flat FRW solutions are related by a "scale-factor dual- 
ity" 0] 

a — > —a , — > (j) — 6a . (3-H) 



This scale-factor duality relates expanding cosmologies to contracting ones. Lidsey []3TJ] has 
shown that this duality can be extended to all values of uj for vacuum solutions in Brans- 
Dicke gravity, 



'2 + 3cA „ / 1 + u 
a — > 1 a 



A + 3ujJ V4 + 3u;, 



which reduces to Eq. ( p.ll|) when u = —1. Thus the fast and slow solutions defined above 



are interchanged under the above duality transform^ 

(p±, Q±) -»• (Pt, ?t) ( 3 - 13 ) 



1 The fast, (p+, q+), and slow, solutions should not to be confused with the (+) and (— ) branches 

defined in the context of low-energy string cosmology by Brustein and Veneziano |33] . 
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4 Autonomous phase-plane 



In the absence of analytic solutions for general FRW models we seek a qualitative description 
of the evolution of the cosmological scale factor and scalar field with respect to time 
2~4|, E5). We shall now show that we can reduce the equations of motion (12.31) and (I2.4D 



We shall now show that we can reduce the equations of motion (|2.3|) and ( j2.4|) to a 
two-dimensional autonomous phase-plane whose qualitative evolution can be described for 
any values of the Brans-Dicke parameter u or barotropic index 7. 

The trace of the energy-momentum tensor for matter acts like a potential gradient term 



driving the scalar field $ in Eq. (2^4) and so, in analogy with the phase-plane analysis of 
self- interacting scalar fields in general relativity |34|, , we can absorb the matter terms 
on the right-hand side of both the field equations Q2.3| ) and (|2.4|) by introducing a new time 
coordinate, r, such that 



dr 



8vrp 



(3 + 2cj)$ 



dt 



We introduce two new variables |24|] 



x = 



V 



1$' 

2~¥ 

a/_ 

a 



Substitution of these variables in the equations of motion (|2.4|) and ( |2.3| ) yields 

4-37 

2 



X 



2 3(2-7) 
-x xy 



-2(1 -3a 2 ) 



y 



3a 2 



-x + 3xy + 



37 



-y 



3(4-37)a 2 + 37-2 
6a 2 



where for convenience we have defined 



(4.1) 

(4.2) 
(4.3) 

(4.4) 
(4.5) 



a 2 = 



2to + 3 

We will consider only models for which uj > —3/2 and thus a 2 > 0. 
The constraint equation (|2.6|) becomes 



{x + yf 



3$A; 



iirpa 2 a 2 



(4.6) 



(4.7) 



By comparison with the Einstein frame constraint equation ( |2.9| ) we see that (x+y) describes 
the expansion in the Einstein frame, and x 2 the effective energy density of the Brans-Dicke 
field relative to the barotropic fluid. The spatially flat (k = 0) models lie on a separatrix in 
the phase-plane: 



x + y = ± 



VTTx 2 

V3< 



a 



Note that in order for k = models to have a turning point in the scale factor (y = 0) we 
require 3a 2 > 1, i.e., to < 0. 
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The equations ( f4.4|) and ( [4.5|) are two first-order differential equations forming a au- 
tonomous dynamical system in x and y. The simultaneous solution of such equations defines 
a set of trajectories in the (x, y) plane. The graphical representation of these equations form 
a system of 'integral curves' where 

dy = Q(x,y) 

dx P{x,y) 1 ' 

where the functions P(x, y) and Q(x, y) are defined by the right hand sides of ( |4.4| ) and ( |4.5| ) 
respectively. This allows the qualitative evolution of the variables x and y to be determined. 



4.1 Critical points in the finite plane 

The key aspect of a phase plane analysis is to determine the critical points (xi,yi) for the 
system. These are fixed points where the time derivatives of the variables vanish, 

x' = , y' = . (4.10) 

For an integral curve approaching such a point, the critical point represents the asymptotic 
behaviour either at the start or the end of the evolution. 

Firstly we seek critical points at finite values (xj, yi) in the phase plane. Requiring x' = 
implies 

Requiring in addition that y' = generates two possible solutions: 

(a) 3a 2 (x a + y a ) 2 = (l + x 2 a ) (4.12) 
or (b) (2 - 37)2/6 = 2x b (4.13) 



Simultaneous solution of Eq. (|4.11| ) with either Eq. (|4.12| ) or ( [4.13|) gives expressions for the 
critical points: 

[x ^ ±) = ±| "JLM , (2 -7) -^(4-37) \ 

/3(2- 7 ) 2 -a 2 (4-3 7 ) 2 a,/3(2- 7 ) 2 -a 2 (4-37)V 



(^) = ±(^^=) (4.1*) 

Note that from the definitions of the variables x and y in Eqs. (|4.2|) and ( f4.3|) the phase 
plane is symmetric under the 180° rotation x — > —x and y — > —y plus time reversal r — > — r. 
Thus for every critical point Xi + ,yi + we have a critical point In the remainder of 

this section we shall refer only to the critical points with Xi + + yi + > which correspond to 
expanding solutions in the Einstein frame. It should be taken as read that there will be a 
point with = — Xi + and = —yi+ which corresponds to the time reversed solutions. 

Also from the definitions of x and y we see that the fixed points (x«, yi) correspond to 
solutions for $ oc e 2XlT and a oc e VtT with respect to the time coordinate r. The evolution 
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at the fixed points is self-similar as the nature of the solutions is unaffected by a rescaling 
t — ► t + At. Integration of Eq. (|4.1|) then gives 



* - t m = - J(Z + **)*o<& e ( 2 ^ + s m )./ 2 . (416) 

This ensures that we can only have cosmic time t —>■ oo at critical points in the finite phase- 
plane when r — ► oo (for 2a;j + 372/j > 0) and t — > — oo when r — ► — oo (for 2xj + 3jyi < 0). 
Critical points thus correspond to power-law solutions with respect to cosmic time for the 
scale factor a and scalar field $ of the form given in Eq. where 



p. = , ft = . (4.17) 

2xi + 3jyi ' 2xi + 3jyi 



4.1.1 Point (a) 



Critical point (x a ,y a ) always lies on the k = separatrix, as can be seen by comparing 
Eqs. ([4TT2D and (g^ . Substituting Eq. flOD into Eq. (pTTTD we see that point (a) corre- 



sponds to Nariai's power-law solution for k = given in Eqs. ( p.2|) and (|3.3| ). 

From Eq. (|4. 14|) we see that point (a) exists for a < v / 3(2 — 7)/|4 — 37 1 which implies 
that oj > ou a (j), where we define 

It can be seen that (x a ,y a ) tends to infinity as oj — > ^(7). In the limit 7 — > 2 we require 
u; — >■ 00 for (x a ,y a ) to exist at finite values in the phase plane. On the other hand, as 
oj — > —3/2 we find that 7 is restricted to 4/3 for point (a) to exist, as shown in Figure |l[ 
For point (a) to correspond to an expanding universe, a > 0, requires y a > and hence, 



from Eq. Q4.14j ), oj > (^+(7), where we define 



W+ ( 7 ) = — ^- (4.19) 
2-7 



Note also that we have t — > +00 as r — > +00 in Eq. (|4.16|) at point (a) only if u; > c<j_(7), 
where we define 

W-(7) = „ ,~ 4 x • (4-20) 
37(2 - 7) 

Conversely we have t — > — 00 as r — > —00 only for < oj- (7). In the particular case 
a; = ^-(7), point (a) corresponds to a non-singular exponential expansion, a oc exp(Ht), 
where the Hubble rate, H, is a constant. 

Nariai's solution corresponds to an accelerated scale factor, a > 0, if p > 1 or p < 0. For 
p < this requires ^-(7) < oj < oj + (j), in which case point (a) corresponds to a contracting 
universe. For p > 1 we require a; > ^(7) when 7 < 2/3, where we define 

- (2 - 7 K 2 -3 7 ) • < 4 - 21 > 
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Figure 1: Lines dividing regions in parameter space with different cosmological behaviour 
at critical point (a). The dashed line represents u) a (j), defined in Eq. ( f4.18|) . The single dot- 



dashed line represents co*(y), defined in Eq. ( |4.21| ). The double dot-dashed line represents 
u + (y), defined in Eq. (^4.19|) . The dotted line represents as defined in Eq. ( |4.20| ). 



or ^0(7) < uj < ^-(7) when 1 < 7 < 4/3. But only in the former case do we have u > ^+(7) 
and hence an accelerated, expanding universe. See Figure [I]. 

For the particular case 7 = (corresponding to a false vacuum energy density), point (a) 
corresponds to 

a oc r +1/2 , $ oc t 2 (4.22) 
which is the basis of models of extended inflation solution [T^] when u > u;*(0) = 1/2. 

4.1.2 Point (b) 

Critical point (xb,yt,) is a novel power-law solution which corresponds to a curvature scaling 
solution [0, |2§. Substituting Eq. ( ftT5D into Eqs. (Q and ( ^T7|) we find 



a oc t , $ oc t 2 -' 3h . (4.23) 

Thus we have $ oc a 2-37 and so the gravitational effect of the matter remains proportional 
to that of the curvature in the constraint Eq. (|2.6| ) 

d\ 2 p k 
- oc — oc — 
$ a 2 



a J 



(4.24) 
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The universe expands with a oc t, as it would in general relativity when 7 = 2/3 and p oc 1/a 2 
in general relativity. In Brans-Dicke theory this solution exists for all 7 < 2/3. Point (b) 
maybe described as marginal inflation, in the sense that although the spatial curvature does 
not vanish with respect to the gravitational energy density in the constraint equation (|2.6|) . 
it does not grow either, and remains proportional to the gravitational energy density as 
t — > 00. 

Point (b) exists for 7 < 2/3, irrespective of the value of w [23, 123]. In the limit 7-^2/3 



we have Xj, — > and y^ — > 00. Note that 2x fe + Syyb > whenever point (b) exists and thus 
we have t — > +00 as r — > +00 in Eq. ( f4.16| ). 

When us = u)*(y), defined in Eq. ( |4.21| ), point b lies on the k = separatrix and the 
critical points (x a , y a ) and (xf,, y^) coincide. When uj > oj* point (b) lies in the k > region, 
while for to < uj* it lies in the k < region. 

4.2 Stability of critical points in the finite plane 

To determine the stability of these solutions we expand about the critical points 

x = Xi + u (4.25) 
V = Vi + v (4.26) 

and substitute into the equations of motion for x and y, Eqs. ( |4.4| ) and (|4.5|). This yields, 
to first-order in u and v , 

„- = _( 2lj+ 5epi^ M _(^^ B (4 , 7) 

/ 4(1 - 3a 2 ) \ . , . , . . 

v = ( 3q; 2 x i + 3yAu + {3x i + (3y-2)y i )v (4.28) 

We seek the two eigenvalues A± such that 

u' = \±u , u' = \±v , (4.29) 

for the eigenmodes u = K±v, where K± are constants. The eigenvalues are given by the 
roots of the quadratic 

A 2 + B\ + C = (4.30) 

where 

10-97 f A \ 

B % = ^—y l -x i (4.31) 

C, , -6 ( t 1 - ^ - ^ + l) - 2(3 7 - 2)^ - 3(2 - ~ 2) ,| (4,2) 

The critical point is a stable, attractor solution if the real part of both the eigenvalues 
is negative (i.e., Bi > and C, > 0). Linear perturbations about such points are then 
exponentially damped. A saddle point in the phase-plane has one negative and one positive 
real root (Cj < 0). An unstable node has eigenvalues with positive real parts (Bi < and 
Ci > 0) corresponding to growing perturbations. 
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For the critical point (6) we find that > for all values of 7 < 2/3 where (b) exists in 
the finite plane. is positive for uj < uj* defined in Eq. (|4.21| ) and point (b) is then a stable 
attractor [^4|, For uj > uj* we find Cb < and point (6) is then a saddle point |24], ^5 . 



For point (a) we find C a > for c<j > c<j*. We also have B a > in this region and thus 
the point (a) is a stable attractor. For uj < uj* we have C a < and the critical point (a) is a 
saddle point. However point (a) is always the attractor for solutions on the k=0 separatrix 
when it exists. 

For the case uj = uj* the equations are degenerate and the critical points (a) and (b) 
coincident. Further analysis is required to determine the nature of the critical point. It can 
be shown that for k < the critical point remains an attractor, whilst for k > the critical 
point is a saddle point. In the phase plane as a whole the point is a col-node. 



4.3 Critical Points at Infinity 



In order to fully describe the qualitative evolution of the system, we must understand the 
nature of the phase plane at infinity. This corresponds to the regime where the dynamical 
effect of the fluid density becomes negligible either with respect to the expansion (y — > ±00) 
or the evolution of the Brans-Dicke field (x — > ±00). 

Compactification can be achieved by projecting the infinite (x, y) plane onto the unit 
disc (r, 6) where 

R 

y 



X 



R cos 9 , y = R sin 9 , 
The equations of motion for r and 9 then reduce to 



1 + R 



r9' = (1 — r) [cos9y' — sin 9 x'\ 
r' = (1 — r) 2 [cos 9 x' + sin 9 y'} 



where x' and y' are given by Eqs. ( (4.4|) and 
At infinity (the limit r — > 1) we find 



(4.33) 



(4.34) 
(4.35) 



9' 



3a 2 (l 



cos 9 



cos 9 



3a 2 (sin 9 + cos 1 



2 (sin 9 + cos 1 



• cos 2 9 
2 



3a 2 



cos 2 9 



3(2 - 7 ) 



sin 9 



(4.36) 
(4.37) 



The equation for 9' determines the behaviour of the solutions at infinity (r = 1) and is 
independent of 7. Solving for 9' = as r — > 1 we find the critical values 9 = 9{, where 



tan 9 C 



-1 + 



V3, 



tan 9 a 



1 



a. 



cos# P = . 







(4.38) 



Clearly there are two solutions for each case in the complete interval (0,2tt). For each 
solution 9 i+ there is a solution 9i_ which is related by 9i_ = 9 i+ + n and corresponds to the 
time reverse of the original solution. Just as we did for the critical points at finite values of 
x and y, in the rest of this section we will consider only the solutions 9 i+ which lie in the 
interval (— 7r/4, +37r/4) corresponding to Xj + yi > 0. 
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Note that \d9/dr\ — > oo whenever 9' 7^ as r — > 1 and thus all trajectories (except those 
with r = 1 at all times) approach the critical points at infinity along a radial trajectory. 
Approaching infinity with 9 = 9i at a finite value of r — > r* we have, from Eqs. (|4.4|) , ( |4.5| ) 
and flOTp , 

$ oc |r - ^l-^Wi) , a oc |r - ^-^/M) . (4.39) 
Integration of Eq. ( j4.1| ) then gives 

t-Uoc\T-T*r m)/m) , (4.40) 

where 

2cos0j + 37 sin ^ -2f(9i) 



F(9 t 



2 COS 9; + 



2 

3-2 (sin 9i + cos 9A 2 + — - cos 2 6{ 

6a z 



sin 0* . (4.41) 



Thus the critical points at infinity also correspond to power-law evolution of the scale factor 
and scalar field with respect to the proper cosmic time t of the form given in Eqs. fl4.17|) 
with 

sin^ 2cos^ 
P F{9 i ) ' q F{9 i ) ■ ^ 



We find t -> +00 as r -> r» in Eq. ( ^40] ) for F(^) > and /(^) > or t -> -00 as 
r -> r, for F(^) < and /(0 4 ) < 0. 

4.3.1 Points (c) and (d) 

Points (c) and (d) lie at either end of the k=0 separatrix, defined by Eq. ( |4.8|) . Substitution 
of tan 9i = — 1 ± 1/ y/3a into Eq. ( 4.42j ) gives the exponents given in Eq. (|3.9| ) for the power- 



law evolution of a and $ with respect to the proper time t. These are spatially flat (k = 0) 
vacuum solutions, first given by O'Hanlon and Tupper p8 |. In terms of the analysis by 



Gurevich et al |17| , points (c) and (d) correspond to the slow and fast solutions respectively. 

Point (d) is always an expanding solution (sin 9 > 0), but point (c) corresponds to a 
contracting solution if uj < 0. 

4.3.2 Point (e) 

Point (e) lies on the x-axis as y — * 00 in the k < region. Substituting cos# e = into 
Eq. (|4.42|) gives the evolution of the scale factor and scalar field with repect to proper time 

t, 

$ — > constant , a oc t (4.43) 
Point (e) is therefore the general relativistic curvature-dominated Milne universe [[2~5f| . 
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4.4 Stability of points at infinity 



By inspecting Eq. (|4.36|) for 9' we see that solutions at infinity (r = 1) approach point (e) 
but diverge away from points (c) and (d) independently of the parameters 7 and uj. The 
overall stability of the points depends upon the sign of r' = f{9.j) approaching these points 
from r < 1. 

For point (e) we have cos^e = and so from Eq. ( |4.37| ) we have f(9 e ) = (37 — 2)/2. Thus 
point (e) is a stable attractor for 7 > 2/3. We also find F(9 e ) = 1 in Eq. (|4.40|) and thus 
t — > +00 as t — > t* at point (e) whenever 7 > 2/3 [p5] . 

For points (c) and (d) we can eliminate l/3a 2 from Eq. (|4.37| ) by the subsitution of 
(1 + tan 9t) 2 using Eqn. ( [4.38[ ), and we obtain 



f(0i) = ~ 3(2 2 7) sm^ - cos9 i (4.44) 
We find that whenever the critical point (a) exists in the finite phase plane [u > £^0(7) 



defined in Eq. ( 4.18 )1 we have f(9i) < for both points (c) and (d) so that the critical 
points are unstable nodes. When point (a) no longer exists the behaviour at the critical 
point depends on the value of 7. For 7 < 4/3 point (c) is a saddle point but is an attractor 
for k = models, whilst (d) is an unstable node. For 7 > 4/3 the converse is true. 

Note also that for points (c) and (d) we find F{9i) = 2 cos 6>j + 3 sin 9i in Eq. ( [4.40|) , we 



have F{9d) > for all uj, whilst F(9 C ) > for uj > —4/3. Thus for point (c) we have t — > +00 
as t — > r* for —4/3 < uj < u a (l) (when 7 < 1), and conversely t — ► —00 as r — » r* for 
co> < —4/3 (when 7 > 4/3) or for LJ a (j) < uj < —4/3 (when 1 < 7 < 4/3). For point (d) 
t — > +00 as r — >• r* for uj < uj a {l) (when 7 > 4/3), there are no regions of the phase plane 
where t — ► —00 as r — > r*. 



5 Qualitative evolution 

In this section we present a qualitative description for the dynamical evolution of the phase 
plane for different regions of the (uj, 7) parameter space. See Figure ^[ 

5.1 Region I: ^(7) < u < 00 , 7 < 2/3. See Figure 0. 

Critical points (a±) and (6±) exist at finite values of x and y, with point (a+) being the stable 
late-time attractor on the x + ?/>0, /c = separatrix. This parameter regime is inflationary 
as point (a + ) corresponds to accelerated expansion and is the attractor solution for all k < 
models with x + y > and a non-zero measure of the k > models. The remaining > 
trajectories and all k < models with 2 + y < approach one of the vacuum solutions (c_) 
or (gL). We find that only those closed (A; > 0) models above a second sepratrix passing 



through point (b + ) reach the inflationary solution at late times [ 24] , |25| . For instance, when 
uj = 18(1 — 7)/(2 — 37) 2 the second separatrix is the straight line y = 2x/ (2 — 7), but we do 
not have an analytic expression for this separatrix in the more general case. 

The generic early-time behaviour is given by the matter dominated solution a_ or the 
vacuum solutions c + or d + , while points b± and e± are saddle points. 
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Figure 2: Parameter regimes defined in Section |5|. The dashed line represents oj a {l) and 
the single dot-dashed line represents a;* (7). 



5.2 Region II: ^(7) < u < a;* (7) , 7 < 2/3. See Figure 

Critical points a± and &± exist at finite values of x and y. We have marginal inflation as 
point b + the stable late-time attractor solution for all k < models with x + y > 0. All 
k > trajectories and /c < models with x + y < collapse to the vacuum solutions c_ or 
d_. 

The generic early-time behaviour is given by the scaling solution 6_ for < models 
with x + y < or the vacuum solutions c + or d + for fc < trajectories with x + y > and 
all k > models. Points a± and e± are saddle points. 

As uj — > u a from above (and we approach region III) point (a + ) moves along the /c = 
separatrix towards point (c + ) at infinity, and (a_) moves off towards On the other 

hand, as 7 — > 2/3 from below (and we approach region IV) we see that points (b±) move off 
towards the curvature dominated solutions at points (e±). 



5.3 Region III: u < cj a (7) , 7 < 2/3. See Figure 0. 

Only points b± now exist at finite values of x and y. We have marginal inflation as point b + 
is the stable late-time attractor for all k < models with x + y > 0. All k > models and 
< trajectories with x + y < approach the vacuum solution gL at late times. Generic 
k — solutions with a; + y > evolve from d + to c + while = solutions with x + y < 
evolve from c_ to gL. 
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e+ 




e- 



Figure 3: Example of phase-plane in region I (7 = 1/3, u — 6). The solid triangle represents 
critical points (a±). The solid square represents the critical points {b±). The dashed line is 
the k = separatrix. 

The generic early-time behaviour is given by the scaling solution b_ for k < models 
with x + y < or the vacuum solution d + for k < trajectories with x + y > and all A; > 
models. Points c± and e± are saddle points. 

As 7 — > 2/3 from below (and we approach region VI) points (6±) tend towards points 

5.4 Region IV: ^(7) < cj < oc , | < 7 < 2. See Figure § and 0. 

Only points a± exist at finite values of x and y. Point e + is the stable late time attractor 
for k < solutions with x + y > 0. All > models and k < trajectories with x + y < 
approach one of the vacuum solutions c_ or gL at late times. Generic k = solutions with 
x + y > evolve from vacuum solutions c + or d + towards the matter dominated solution a + . 

The generic early-time behaviour is given by point e_ for k < models with x + y < 
or the vacuum solutions c + or d + for < trajectories with x + y > and all k > models. 
Points a± are saddle points. 

For 7 < 4/3 we see that as u — > uo a from above (and we approach region V) point (a+) 
moves off towards point (c+), and (a_) moves off towards (cL) at infinity. On the other 
hand, as u; — > uo a from above for 7 > 4/3 (and we approach region VI) point (a + ) moves off 
towards point (d+), and (a_) moves off towards (c_) at infinity. 
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5.5 Region V: u < u a (j) , 2/3 < 7 < 4/3. See Figure g. 

There are no critical points at finite values of x and y. Point e + is the stable late time 
attractor for k < solutions with x + y > 0. All > models and k < trajectories with 
x + y < approach the vacuum solution oL at late times. 

The generic early-time behaviour is given by point e_ for k < models with x + y < or 
the vacuum solutions d+ for k < trajectories with x + y > and all k > models. Points 
c± are saddle points. 

5.6 Region VI: u < Lu a (~f) ,4/3 < 7 < 2. See Figure 0. 

There are no critical points at finite values of x and y. Point e + is the stable late time 
attractor for k < solutions with x + y > 0. All k > models and /c < trajectories with 
x + y < approach the vacuum solution c_ at late times. 

The generic early-time behaviour is given by point e_ for k < models with x + y < 
or the vacuum solutions c + for A; < trajectories with x + y > and all k > models. 

Points d± are saddle points. 

6 Non-singular evolution 

In general relativity it is well-known that all spatially flat or open (k < 0) FRW models with 
a barotropic fluid have a semi-infinite lifetime and posseess a singularity either at a finite 
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cosmic time in the past (for expanding models) or in the future (contracting solutions) so 
long as 7 > 0. Closed models (k > 0) have a finite lifetime with both a past and a future 
singularity unless 7 < 2/3 in which case non-singular evolution, i.e., infinite proper lifetime, 
is possible. 

It is natural to ask whether in the context of Brans-Dicke gravity it is possible to have a 
non-singular universe in open or flat models, and/or for 7 > 2/3. Our phase-plane analysis 
provides a qualitative description of the asymptotic evolution that is possible in the different 
parameter regimes in terms of the time parameter r, and equations ( (4.16| ) and ( |4.40| ) relate 
the asymptotic behaviour of r to the proper cosmic time t. 

As oj — > 00 the critical points (a + ) and (b + ) with Xi + yi > follow the pattern described 
above for expanding models in general relativity with semi-infinite proper lifetimes. Point 
(&+), whenever it exists (7 < 2/3), always describes a solution with a past singularity, but 
continuing in to the infinite future. Point (a + ) similarly describes a semi-infinite lifetime 
solution with past singularity for uj > max{u; a , Points (a_) and (&_) are just the time- 
reversed solutions of points (a+) and (b + ) respectively, and thus are also semi- infinite but 
with the future singularities. This suggests that solutions which interpolate between points 
(a/b-) and (a/b + ) would have an infinite lifetime, but the phase-plane plots show that the 
solutions which cross the line x + y = always start from a critical point with Xi + yi > 
and end up at a critical point Xj + yj < which would yield a finite lifetime. These are 
analogous to the closed models in general relativity which evolve from a past singularity to 
a future singularity. 

However for u a < uj < uj_ (which is possible for 1 < 7 < 4/3) the solution at point (a + ) 
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Figure 6: Example of phase-plane in region IV when 7 < 4.3 (here 7 = l,u — 6). Symbols 
as for Figure ||. 



extends to the infinite past and the singularity lies in the future, whereas it is the point (a_) 
which has a singularity in the past. But in this case we find no trajectories which link point 
(a + ) to point (a_). (Points (b±) do not exist in this parameter regime.) 

Solutions at the fixed points (a) or (b) are necessarily singular to the past or future with 
semi-infinite lifetime because they correspond to power-law evolution for the scale factor, 
except when uj = uJ a {l)- In this particular case the scale factor expands exponentially and 
the solutions at points (a±) is indeed non-singular. However in this case the solution is 
unstable as points (a±) are saddle points in the phase-plane. 

Generic solutions in the phase-plane interpolate between two of the critical points includ- 
ing at least one of the vacuum solutions (c±), (d±) or (e±) at infinity. As uj — ► 00 points (c+) 
and (d + ) are early-time attractors [fO < in equation (|4.37|) 1 as r — > r* and this corresponds 
to a singularity at finite time t* [F(9) > in equation (|4.40| )1. For 7 < 2/3 point (e+) is 
also an early-time attractor with a singularity at a finite time £*, but for 7 > 2/3 point (e + ) 
becomes a late-time attractor where t — > +00 as r — ► r*. For all values of 7 we are unable 
to construct non-singular trajectories as uj —>■ 00. The situation is actually worse than in 
general relativity plus barotropic fluid due to the presence of the scalar field. Even in the 
limit uj — * 00 the Brans-Dicke field leads to singular evolution at early or late times even in 
closed models when 7 < 2/3. 

However for ^(7) < uj < —4/3, which is possible for 7 > 1, we have t — ► — 00 at 
point (c + ) which is still an early-time attractor. Thus trajectories in this parameter regime 
that originate at point (c + ), which has no past singularity, are non-singular if they connect 
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Figure 7: Example of phase-plane in region IV when 7 > 4.3 (here 7 = 1.5, u — 6). Symbols 
as for Figure ||. 



to a critical point with no future singularity. We find three parameter regions in which 
non-singular trajectories occur: 

1. UJ a < UJ < 0J_ 

• Open models (k < 0) can originate at point (c + ) in the infinite past and approach 
point (e + ) in the infinite future. 

2. max{w_, u a } < uj < —4/3 

• Open models (k < 0) can originate at point (c + ) in the infinite past and approach 
point (e + ) in the infinite future. 

• Flat models (k = 0) can originate at point (c + ) in the infinite past and approach 
point (a + ) in the infinite future. 

• Closed models (k > 0), with 7 > 4/3, can originate at point (c + ) in the infinite 
past and approach point (c_) in the infinite future. 

3. uj < min{c<j_, — 4/3} 

• Open models (k < 0) originate at point (c + ) in the infinite past and approach 
point (e + ) in the infinite future. 
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• Flat models (k = 0) originate at point (c+) in the infinite past and approach point 
(d + ) in the infinite future. 

• Closed models originate at point (c+) in the infinite past and approach point (c_) 
in the infinite future. 

Note that the time reverse of these solutions linking points (e_) and (c_), etc., will also be 
non-singular in the appropriate parameter regime. 

We should stress that an infinite proper lifetime is a necessary, but not sufficient, condition 
for these solutions to be non-singular. In particular, Kaloper and Olive [^] have recently 
emphasised that as the metric is, by definition, only minimally coupled to the other fields in 
the Einstein frame, gravitational waves follow geodesies in this frame. Thus if the solutions 
have a singularity in the finite past in terms of the proper time in the Einstein frame, defined 
in Eq. ( |2.8|) , they may still be geodesically incomplete. One can show that all the critical 
points with Xi + jji > have singularities at a finite proper time in the past, and extend into 
the infinite future, in the Einstein frame. 

7 Summary 

We have constructed an autonomous phase-plane describing the evolution of the scale factor 
and Brans-Dicke field for FRW cosmologies containing a barotropic fluid in Brans-Dicke 



20 




gravity. We have improved upon previous phase-plane analyses |25| by presenting all the 
FRW models in a single phase-plane. We are thus able to present a qualitative analysis 
of the general evolution of homogeneous and isotropic cosmologies in Brans-Dicke gravity, 
recovering previously known power-law solutions as fixed points in the phase-plane. These 
fixed points correspond to self-similar evolution [23|. This is possible due to the scale- 
invariance of the Brans-Dicke gravity theory which contains no characteristic value for the 
Brans-Dicke field, $. 

We find four critical points at finite values of the phase-plane when 7 < 2/3 and u > ^(7) 
[see Eq. ([4.18 )]. Two of the fixed points correspond to Nariai's power-law solutions [p9| |, 
a oc t p , $ oc t q , for spatially flat models (k = 0). One of these fixed points is expanding 
y > and the other is contracting y < 0. The remaining two fixed points describe a novel 
self-similar curvature-scaling evolution for spatially curved models with a oc t and 

$oct 2 - 3 ^. 

Nariai's expanding solution is only an attractor at late times if uj > ^(7) [see Eq. ( |4.21| )]. 
This corresponds to power-law inflationary solutions with exponent p > 1 and a non-zero 
measure of spatially curved models approach this flat-space solution. For instance, for a false 
vacuum energy density (7 = 0) we have uj* = 1/2. On the other hand, for a given value of 
uj the necessary condition for inflation becomes 7 < 7*(cu), where [0, 24], from Eq. ( [1.21| ), 



7* [uj ) 




(7.1) 
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which is always stronger than in general relativity (u — > oo) where 7* = 2/3. For 7*(cj) < 
7 < 2/3 the expanding curvature-scaling solution with k < is a late-time attractor for 
k < models. 

For 7 > 2/3 and u > u a {j) only the two critical points corresponding to Nariai's solutions 
remain at finite values in the phase-plane, while for 7 < 2/3 but u> < u a (j) only the curved 
space fixed points remain. For 7 > 2/3 and uo < oo a {l) their are no fixed points in the finite 
phase-plane. 

By projected the infinite phase-plane onto a unit disc we studied the asymptotic behaviour 
of solutions at infinity where the fluid density becomes negligible. Here we recover vacuum 
solutions corresponding to the slow and fast branches of O'Hanlon and Tupper's vacuum 
solutions J28| for k = 0, and the Milne universe for k < models. The expanding Milne 



universe is a late-time attractor (and the collapsing solution is an early-time atractor) for 
7 >2/3. 

We have seen that it is only possible to construct solutions with infinite proper lifetimes 
for values of the Brans-Dicke parameter uj < —4/3. However, in all cases the solutions 
possess a curvature singularity in the conformally related Eisntein frame. 
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